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摘 要 : LAE KEK MAB RL (SPED) 作为 评估 指标 ,基于 渭河 流域 28 个 气象 站 点 1961—2017 年 
实测 降水 量 和 气温 数据 ,采用 Mann-Kendall( M-K) 趋势 检验 、 经 验 正 交 函 数 以 及 小 波 变换 等 方法 分 
析 渭 河流 域 干旱 时 空 变化 特征 ,并 研究 渭河 流域 干旱 与 6 种 大 尺度 气候 因子 之 间 的 相关 关系 , 进 一 
步 探 讨 主要 气候 因子 对 流域 干旱 时 空 分 布 特征 的 潜在 影响 。 研 究 表 明 : 渭 河流 域 在 1961 一 2017 年 
间 整 体 呈现 出 变 旱 的 趋势 。 通 过 经 验 正 交 耳 数 分 解 ,渭河 流域 干旱 分 布 场 主要 有 3 种 典型 模 态 类 


Al ,分别 为 全 局 型 


\ 西 北 一 东南 反 向 分 布 型 以 及 东 一 西 反 向 分 布 型 ;同时 ,大 尺度 气候 因子 南方 涛 动 


指数 SO] 与 流域 干旱 分 布 场 具 有 更 好 的 相关 关系 ,对 该 区 域内 干旱 变化 有 较 强 的 影响 。 
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受 全 球 气候 变 暧 影响 ,干旱 已 成 为 全 球 范围 内 
最 复杂 、 最 常见 的 自然 灾害 之 一 ， 。 特 别 是 在 水 资 
源 紧缺 的 地 区 ,频繁 的 干旱 已 严重 降低 了 人 们 的 生 
产生 活水 平 ,甚至 威胁 到 人 类 的 财产 及 生命 安全 ,对 
社会 的 发 展 带 来 不 利 影响 。 在 此 情况 下 ,开展 对 干 
旱 时 空 分 布 特征 的 科学 研究 ,以 达到 防御 旱灾 、 降 低 
干旱 影响 的 目的 ,已 经 成 为 国内 外 研究 的 热点 ” 。 
我 国 北 方 地 区 水 资源 紧缺 ,同时 受 季风 气候 的 影响 ， 
北方 地 区 发 生 干 旱 事 件 较 为 频繁 ”。 渭 河流 域 位 


渭河 流域 ; 干旱 ; SPE; 经 验 正 交 函数 ; 小 波 变换 


散 指数 (Standardized Precipitation Evapotranspiration 
Index ,SPEI) , SPEI ERLA T PDSI 对 蒸 散 的 响应 ， 
也 融合 了 SPI 空间 一 致 性 多 时 间 尺 度 的 优点 呈 , 因 
此 , 随 着 SPEI 的 提出 ,许多 学 者 使 用 该 指数 对 不 同 
区 域 做 了 不 少 研究 ”。 例 如 , 李 忆 平等 研究 指出 
SPEI 指数 考虑 了 温度 的 因素 ,对 全 球 气温 快速 上 升 
导致 干旱 化 反映 更 加 敏感 ; 张 玉 静 等 "" 应 用 SPEI 
对 华北 地 区 干旱 特征 进行 研究 ,表明 SPET 对 干旱 的 
表征 准确 ,适用 于 对 流域 干旱 的 评估 与 检测 ; 张 岳 军 


于 我 国 半 湿 润 半 干旱 过 渡 的 北部 地 区 ,干旱 灾害 频 
繁 ,流域 内 的 社会 经 济 发 展 已 经 受到 干旱 的 严重 影 
响 中 ,因此 ,进行 关于 渭河 流域 干旱 的 研究 具有 重 

目前 ,国内 外 对 不 同 干旱 程度 的 量化 已 相继 提 
出 多 种 干旱 指标 "" ,其 中 较为 成 熟 的 主要 有 帕 默 尔 
干旱 指数 (51 (Palmer Drought Severity Index, PDS/) 
和 标准 化 降水 指数 "1 ( Standardized Precipitation In- 
dex, SPI), 2010 年 ,VICENTE-SERRANO 等 5 在 标 
准 化 降水 指数 SPI 的 基础 上 ,增加 考虑 蒸 散 作 用 对 
干旱 的 影响 ,提出 一 种 新 的 干旱 指标 ,标准 化 降水 落 


O 收 稿 日 期 : 2019-05-24; 修订 日 期 : 2019 -08 -27 


等 '” 通过 研究 表明 SPEL 可 以 更 灵活 地 反映 月 季 干 
旱 变化 特征 , 且 SPE7 与 部 分 大 尺度 气候 因子 存在 较 
好 的 相关 关系 。 

鉴于 此 ,本文 以 标准 化 降水 蒸 散 指数 SPET 作为 
指标 ,采用 Mann-Kendall ( M-K ) 趋势 检验 .经 验 正 交 
函数 ( Empirical Orthogonal Function , EOF) 和 小 波 变 
换 等 分 析 方法 ,重点 从 分 布 场 的 角度 分 析 研 究 渭河 
流域 干旱 的 时 空 变化 特征 ,揭示 其 变化 规律 ;同时 ， 
考虑 各 种 大 尺度 气候 因子 的 可 能 影响 ,识别 影响 渭 
河流 域 干 旱 变化 的 主要 气候 因子 ,并 剖析 主要 因子 
对 流域 干旱 时 空 变化 特征 的 潜在 影响 ,以 期 为 渭河 
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流域 干旱 
论 依据 。 


1 研究 区 概况 及 数据 来 源 


1.1 研究 区 概况 

渭河 流域 发 源 于 甘肃 省 定西 市 渭 源 县 ,于 陕西 
JA PATH WHEEL A E ,地 理 位 置 在 103"57' ~ 
110°16’ E, 33°42’ ~ 37°24'N 之 间 , 总 面积 约 为 
13.48 x 10° km 。 渭 河流 域 总 体 上 呈现 西高 东 低 的 
态势 ,西部 地 区 以 黄土 拷 陵 沟 蜜 区 为 主 ,而 东部 地 区 
主要 是 关中 平原 区 。 流 域 主要 呈 大 陆 性 气候 特征 ， 
夏季 高 温 多 雨 ,冬季 寒冷 干燥 ,年 均 气 温 6~14 C, 
年 均 降 水 量 450 ~700 mm”!  。 滑 河 支流 较 多 ,左岸 


灾害 防治 和 水 资源 综合 利用 提供 重要 的 理 


Administration, NOAA) 地球 系统 人 研究 3 
据 (https://www. esrl. noaa. gov/ ) 。 
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2 研究 方法 


2.1 标准 化 降水 蒸 散 指数 

标准 化 降水 燕 散 指数 (SPE7) 是 VICENTE-SER- 
RANO 等 中 在 标准 化 降水 指数 (SP1) 基础 上 考虑 蒸 
散 作用 而 提出 的 气象 干旱 指标 。SPEI 具有 SPI 空 
间 一 致 性 和 多 时 间 尺 度 优 点 的 同时 ,采用 降水 量 与 
潜在 蒸 散 量 的 差 值 作 为 水 分 慌 亏 代替 SPL 计算 中 的 
降水 量 , 通 过 考虑 地 表 莹 散 对 干旱 的 影响 ,能够 对 全 
球 气温 上 升 导致 的 干旱 化 反应 更 加 敏感 。 因 此 ， 
随 着 SPEI 的 提出 ,该 指数 在 国内 外 研究 中 得 到 广泛 


的 北 洛 河和 泾 河 是 最 大 的 2 条 支流 ,控制 流域 面积 
分 别 约 为 2.70 km? 和 4.54 km? ,分 别 占 渭河 流域 总 
面积 的 20% All 33.7% 。 
1.2 数据 来 源 

本 文 研究 主要 使 用 到 渭河 流域 气象 站 点 实测 数 
据 以 及 全 球 大 尺度 气候 因子 数据 。 气 象 站 点 数据 来 
自 于 我 国 气象 科学 数据 服务 网 , 共 包 括 渭河 流域 及 
其 周边 28 个 气象 站 点 (图 1)1961 一 2017 年 逐 月 降 
水 量 和 气温 数据 ;大 尺度 气候 因子 主要 使 用 到 表征 
ENSO 事件 的 南方 涛 动 指数 (SO7) 太平 洋 年 代 际 振 
荡 (PDO) 、 北 太平洋 遥 相 关 指数 (NP) .北大 西洋 振 
H (NAO) 北极 涛 动 (40)、 大 西洋 年 际 振荡 指数 
(AMO )6 个 主要 类 型 ,数据 均 来 自 于 美国 国家 海洋 
和 大 气管 理 局 (the National Oceanic and Atmospheric 
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图 1 渭河 流域 及 气象 站 点 分 布 示意 图 
Fig.1 Topographic map of the Weihe River Basin 


and meteorological stations 


应 用 。 本 文选 择 以 同时 考虑 降水 量 和 气温 的 SPEI 
作为 指标 , 进行 渭河 流域 干 湿 时 空 变化 特征 分 析 
研究 。 

SPEI 对 干旱 过 程 的 研究 有 不 同时 间 尺 度 ,本 文 
以 12 个 月 为 时 间 尺 度 ,通过 计算 渭河 流域 各 气象 站 
点 SPEI 序列 ,研究 渭河 流域 年 际 干 湿 时 空 变化 特 
征 。SPEI 干旱 指数 具体 计算 过 程 可 见 参考 文献 上 ， 
其 中 对 潜在 蒸 散 的 计算 使 用 以 气温 为 输入 的 桑 斯 威 
特 (Thornthwaite ) 公式 [16] 。 


SPEI 干旱 分 级 情况 如 下 表 1 Stas” 。 


表 1 SPE 干旱 等 级 划分 标准 
Tab.1 Drought classification standard of SPEI 


等 级 类 型 SPEI 
1 无 旱 SPEI > -0.5 
轻 旱 -1.0 < SPEIs -0.5 
3 pE -1.5 < SPEIS -1.0 
4 重 旱 -2.0<SPEIs -1.5 
5 特 旱 SPEI< -2.0 


2.2 Mann-Kendall( M-K ) 趋势 检验 

研究 采用 曼 肯 德尔 秩 和 检验 ( Mann-Kendall 
test, M-K) 进行 时 间 序 列 的 趋势 变化 分 析 '"”"。 在 使 
用 M-K 过 程 中 ,如 果 原 始 待 检 时 间 序 列 自 相关 性 较 
高 ,M-K 检验 结果 往往 会 存在 一 定 误差 。 因 此 ,在 
使 用 M-K 方法 进行 趋势 检验 之 前 ,需要 先 剔除 原始 
时 间 序 列 的 自 相关 性 。 本 文采 用 预 置 白 方法 
(the trend-free pre-whitening, TFPW ) 对 时 间 序 列 进 
行 预 处 理 ,以 剔除 原始 序列 的 自 相 关 性 ,TFPW 预 处 
理 方式 如 下 。 
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对 于 原始 时 间 序 列 |x;,i=1,2,…,n| ,首先 计 
算 其 一 阶 自 回归 系数 : 


_ Cov(%; ,Xi+1 ) 
n Var(x;) 
n—-2 i+l —%) 
= (1) 
EE 3 
采用 以 下 公式 剔除 原 序 列 的 自 相 关 性 : 
X; =X; Y1; (2) 
用 新 的 序列 |x,i=1,2,…,n| 代替 原始 序列 
i ee 


2.3 经 验 正 交 函 数 分 析 法 

经 验 正 交 函数 分 析 ( Empirical Orthogonal Func- 
tion, EOF) 同 时 也 称 为 主 成 分 分 析 (Principal Com- 
ponent Analysis, PCA) ,通过 对 要 素 场 进行 时 空 分 
解 , 提 取 和 矩阵 数据 的 主要 特征 向 量 , 可 以 实现 对 原始 
数据 空间 和 时 间 变 化 特征 进行 综合 分 析 。EOF 分 
析 主 要 原理 与 算法 过 程 如 下 。 

(1) 首先 ,对 空 = 间 数 据 做 预 处 理 ， 得 到 一 站 点 个 
数 为 m\ 时 间 系 列 长 度 为 n BREE X nxn o TE FEE X 
及 其 转 置 X Wy AR 


C, olxxx? (3) 


mxm n 
(2) 计算 方 阵 Cx 的 特征 根 (A ,…， 
征 向 量 Vazn 
cae x am = Voisin X E mxm (4) 


式 中 :矩阵 Vn xm Fe E HE LE, EA BI A 5 UE EAS 
PIEI E 3E mm AY me xm FEAT FAE RE 


Àn) 与 特 


A, 0 0 
0 A, … 0 
帮 = (À ZA SAn Z0) (5) 
0 0 WEN Àn 
SRA A, Aas ,A 分 别 为 方 阵 Co, 的 特征 根 。 
(3) 对 原始 和 矩阵 工 做 投影 ,计算 各 特征 向 量 对 
应 的 时 间 系 数 PC, (也 称 为 主 成 分 ) : 
PG x = Vy XX mxn (6) 


JERE X IKN: 


X xn = axi x PC nxn (7) 


式 中 :PC 为 时 间 系数 矩阵 ; Vi cm EERE ,的 转 
HEE, 

(4) 2—4 , rR RIE m E AY 77 RE AE 
阵 X 的 方差 大 小 用 特征 根 大 小 表示 ,特征 根 越 大 ， 
对 应 模 态 对 空间 要 素 场 的 作用 也 越 大 。 第 有 个 模 态 
的 方差 贡献 率 表示 为 : 


À; 
+ — x 100% (8) 


i=1 


SEPA, A, 分 别 为 排序 为 bi 的 特征 根 ;m 为 站 点 
个 数 。 

(5) 最 后 ,对 最 终 空间 模 态 进行 随机 性 检验 分 
析 。 根 据 相关 研究 ”1 ,在 95% 置信 度 水 平 下 , 特 


征 根 的 误差 为 : 


2 
AA, =A [ae (9) 

式 中 :NM 是 有 效 自由 度 。 如 果 临 近 2 个 特征 的 误差 
范围 未 重 登 , 则 认为 对 应 模 态 通过 检验 ,满足 随机 性 
的 要 求 ;反之 , 则 表明 对 应 空间 模 态 未 通过 显著 性 检 
验 。 
2.4 小波 变换 分 析 
2.4.1 连续 小 波 变 换 ”连续 小 波 变 换 ( Continuous 
Wavelet Transform，CWT) 可 以 对 研究 时 间 序 列 在 时 
域 上 存在 的 多 层次 时 间 结 构 和 局 部 特征 进行 提取 ， 

分 辩 非 平稳 时 间 序 列 的 时 频 局 部 变化 情 
况 -”。 连 续 小 波 变 化 类 型 较 多 ,本 文采 用 在 水 文 
气象 中 运用 较 多 的 Morlet 小 波 , 分 析 渭 河流 域 SPEI 
时 间 序 列 及 EOF 分 析 时 间 系 数 的 周期 特性 。Morlet 
小 波 函 数 形式 为 : 


W(t) sa ete Pe (10) 
式 中 :i 为 时 间 ;w。 是 无 量 纲 频率 ,一般 取 wo =6。 
对 于 时 间 序 列 了 = (xi yay, ) ,连续 小 波 变 
换 定义 如 下 : 


Ws) = /Ey 


RP: Wi (s) 是 小 波 变换 系数 ;6t 是 时 间 步 长 ;s 是 比 
例 因 子 ;n 为 时 间 系 列 长 度 。 

2.4.2 交 又 小 波 变 换 ”交叉 小 波 变换 ( Cross-wave- 
let Transform, XWT) 同时 结合 了 小 波 变换 和 与 交叉 
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谱 分 析 , 可 以 研究 2 个 非 平 稳 时 间 序 列 之 间 的 时 频 
共振 周期 与 位 相关 系 5-s] 。 本 文 主要 使 用 到 交叉 
小 波 变 换 分 析 研究 渭河 流域 干旱 变化 特征 对 大 气 环 
流 指数 的 响应 特征 。 

EWE (s) A WY Cs) SPH SEI TAL FEI X = (x， 
iy 507 ) AY = (yy Vn ，…,y) 的 小 波 变换 系数 , 则 
和 了 之 间 交 叉 小 波 变 换 定 义 为 3: 

WE (s) =Wi(s) x WY’ (s) (12) 

Ap: W (s) 表 示 W (s) HAHH. 
3 结果 及 分 析 
3.1 标准 化 降水 蒸 散 指数 的 时 空 分 布 特征 


基于 渭河 流域 气象 站 点 逐 月 降水 量 和 气温 数 
据 , 以 12 月 为 时 间 尺 度 ,分 别 计算 各 站 点 表征 年 际 


CC SPEI 


SPEI 


5 a 滑动 平均 


干旱 变化 的 SPELT 时 间 序 列 ,并 用 泰 森 多 边 形 法 推 求 
整个 渭河 流域 SPEIL 指标 值 ,图 2 和 图 3 分 别 显示 了 
渭河 流域 SPEI 在 时 间 和 空间 上 的 趋势 变化 特征 。 
如 图 2, 渭河 流域 SPEI 整体 呈 下 降 趋 势 , 上 且 通 过 置 
信 水 平 0.05 显著 性 检验 ( < - 1.96), FEE, SPEI 
平均 值 从 1961 一 1990 年 的 0.34 降低 到 1991 一 2017 
年 的 -0.37; 根 据 SPEIL 旱 漠 分 级 标准 ,流域 最 干旱 
年 发 生 在 1997 年 , SPEI 值 为 -2.05, 属 于 特 早 情 
形 。 如 图 3 ,流域 绝 大 多 数 气 象 站 点 SPEI 时 间 序 列 
值 均 呈 下 降 趋 势 , 且 其 中 13 个 站 点 呈 显著 性 下 降 ， 
可 以 说 明 在 渭河 流域 整个 空间 分 布 上 均 具有 变 早 的 
趋势 。 
3.2 渭河 流域 干旱 时 空 变化 特征 

应 用 EOF 方法 对 渭河 流域 SPEIL 要 素 场 作 时 空 
分 解 , 并 进行 North 显 著 性 检验 一 ,得 到 流域 SPEJ 
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图 2 渭河 流域 SPET 时 间 序 列 
Fig.2 Time series of SPEI in Weihe River Basin 
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图 3 渭河 流域 SPEL 空间 趋势 变化 特征 
Fig.3 Spatial trend characteristic of SPEI in 
Weihe River Basin 
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主要 空间 分 布 模 态 。 如 表 2 , 渭河 流域 干旱 空间 分 
布 前 5 个 空间 模 态 的 特征 值 方差 累计 贡献 率 为 
79.24% ;其 中 ,前 3 个 特征 根 累 计 贡 献 率 达 
72.85% , 且 特征 值 误差 范围 不 产生 重 二 ,通过 North 


表 2 SPEI 序列 EOF 分 解 模 态 的 贡献 率 
Tab.2 Contribution of the EOF mode of SPEI 


in Weihe River Basin 


模 态 特征 什 方差 贡献 Rib He 特征 值 特征 值 
率 /% 贡献 率 / % 误差 下 限 ”误差 上 限 
1 14.16 50.43 50.43 10.37 17. 94 
2 4.87 17.35 67.79 3.57 6.17 
3 1.42 5.06 72.85 1.04 1.80 
4 0.96 3.41 76.26 0.70 1.21 
5 0.83 2.97 79.24 0.61 1.06 


显著 性 检验 ,满足 随机 性 要 求 。 因 此 ,可 以 认为 前 3 
个 空间 模 态 可 用 于 表示 渭河 流域 干旱 的 主要 时 空 变 
化 特征 。 

BRAS 1 特征 向 量 的 方差 贡献 率 为 50. 43% , 远 
高 于 其 他 模 态 ,说明 流域 干旱 的 空间 分 布 特征 主要 
受 模 态 1 的 影响 。 模 态 1 时 空 变 化 特征 以 及 第 一 主 
成 分 (PC1 ) 存 在 的 周期 特点 如 图 4 所 示 。 如 图 4a， 
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图 4 EOF1 空间 分 布 和 时 间 系 数 PCL 及 PCL 连续 小 波 变换 功率 谱 


Fig.4 Spatial pattern of the EOF1, time series of PC1 and continuous wavelet power spectrum of PC1 


在 整个 渭河 流域 , 模 态 1 的 特征 值 均 为 正 值 , 表 明 在 
1961 一 2017 年 渭河 流域 的 干 湿 变化 具有 高 度 的 一 
致 性 , 即 在 模 态 1 影响 下 ,渭河 流域 要 么 整体 干旱 ， 


南 地 区 湿润 ;或 者 西北 地 区 湿润 ,而 东南 地 区 干旱 。 
如 图 Sb, BEAS 2 时 间 系 数 的 趋势 斜率 为 正 ,说 明 模 
态 2 影响 下 渭河 流域 西北 地 区 呈 湿 润 的 趋势 ,而 东 


要 么 整体 湿润 。 此 外 ,特征 值 高 值 中 心 主要 在 渭河 
流域 中 部 区 域 ,说 明 该 地 区 的 干 湿 变 化 程度 更 大 , 且 
在 全 流域 干旱 时 ,中 部 地 区 干旱 程度 更 为 显著 。 如 
图 4b 模 态 1 对 应 时 间 系 数 PC1 显示 ,在 模 态 1 下 ， 
渭河 流域 整体 上 呈现 变 早 趋势 ,其 中 最 湿润 年 份 发 
生 在 1964 年 ,最 干旱 年 份 发 生 在 1997 年 。 通 过 对 
时 间 系 数 PC1 的 连续 小 波 变 换 分 析 , 在 95% 置 信 水 
平 下 PCL 的 周期 特征 并 不 明显 , 仅 在 1965—1972 
年 间 有 2 a 左右 周期 ,以 及 在 1982 一 1990 年 间 有 4 a 
左右 的 周期 (图 4c) 。 

BRAS 2 特征 向 量 的 方差 贡献 率 为 17. 35% , 同 
样 是 影响 渭河 流域 干旱 空间 分 布 的 主要 模 态 之 一 ， 
其 时 空 变化 特征 及 第 二 成 分 (PC2 ) 存 在 的 周期 特性 
如 图 5 所 示 。 如 图 5a, 渭 河流 域 模 态 2 的 特征 值 有 
正 数 也 有 负数 ,在 西北 地 区 为 正 值 , 在 东南 地 区 呈现 
为 负 值 。 在 模 态 2 影响 下 ,渭河 流域 干 湿 呈现 西 
北 一 东南 反 向 分 布 的 态势 , 即 要 么 西北 地 区 干旱 , 东 


南 地 区 有 更 加 干旱 的 趋势 。 同 时 ,图 5b 时 间 系 数 
PC2 的 连续 小 波 变换 结果 显示 ,在 95% 置信 水 平 
下 ,PC2 在 1970—1990 年 间 具 有 时 间 尺 度 4~10 a 
的 周期 ,以 及 在 2000—2010 年 间 有 时 间 尺 度 4~Sa 
的 周期 特性 。 

模 态 3 特征 向 量 的 方差 贡献 率 为 5.06% ,其 时 
空 变化 特征 及 第 三 主 成 分 (PC3 ) 存在 的 周期 特点 如 
图 6 所 示 。 如 图 6a, 模 态 3 的 特征 向 量 在 渭河 流域 
东部 地 区 为 正 值 ,而 在 西部 地 区 为 负 值 , 在 模 态 3 影 
响 下 ,渭河 流域 干 湿 呈 东 一 西 反 向 分 布 , 即 要 么 东部 
地 区 湿润 ,西部 地 区 干旱 ;或 者 东部 地 区 干旱 、 西 部 
地 区 温润。 如 图 6b , 模 态 3 时 间 系 数 具 有 降低 的 趋 
势 ,说 明 在 模 态 3 下 渭河 流域 具有 东部 地 区 逐渐 变 
旱 而 西部 地 区 变 得 湿润 的 趋势 。 图 6c 时 间 系 数 
PC3 的 连续 小 波 变换 ,在 95% 置信 水 平 下 , PC3 在 
1965 一 1972 年 和 2004 一 2012 年 具有 时 间 尺 度 2 ~ 
4a 的 周期 , 以 及 在 1995 一 2010 年 间 有 约 6~10 a 的 
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图 5 EOF2 空间 分 布 和 时 间 系 数 PC2 及 PC2 连续 小 波 变换 功率 谱 
Fig.5 Spatial pattern of the EOF2, time series of PC2 and continuous wavelet power spectrum of PC2 
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(c) PC3 连 续 小 波 变换 功率 谱 
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图 6 EOF3 空间 分 布 和 时 间 系 数 PC3 及 PC3 连续 小 波 变换 功率 谱 


Fig.6 Spatial pattern of the EOF3, time series of PC3 and continuous wavelet power spectrum of PC3 
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分 看 等 :基于 SPE 的 渭河 流域 干旱 时 空 变化 特征 分 析 


周期 特性 。 
3.3 干旱 指数 对 大 尺度 气候 因子 的 响应 
计算 得 到 渭河 流域 年 时 间 尺 度 SPEI 及 EOF 分 
析 时 间 系 数 PCs( PC1、PC2 和 PC3 ) 与 大 尺度 气候 
因子 (SO PDO NP NAO 40 Fil AMO ) 之 间 的 相关 
系数 ,如 表 3 所 示 。 研 究 采 用 上 统计 量 检验 相关 性 
显 车 水平, 时间 序列 长 度 一 共 57 a, 在 显著 性 水 平 
5% 下 ,临界 相关 系数 约 为 0.255。 表 3 RH , SPEI 
及 PCs 与 表征 ENSO 事 件 的 南方 涛 动 指数 S01 的 相 


表 3 气候 因子 与 SPEI 及 PCs 相关 系数 
Tab.3 Correlation coefficients between the climate 
indices and SPEI ,PCs 
类 型 SOI PDO NP NAO AO AMO 


SPEI 0.49 -0.17 0.11 -0.16 0.05 -0.24 
PCI 0.41 0.16 -0.12 0.16 -0.06 0.26 
PC2 0.37 0.05 0.05 0.06 0.17 -0.03 
PC3 -0.28 -0.07 -0.09 0.16 0.09 0.05 


WŒ SOL 为 南方 涛 动 指数 ;PDO 为 太平 洋 年 代 际 振荡 ;NP 为 北 太 平 洋 


HEAD CHAR; NAO 为 北大 西洋 振荡 ;40 为 北极 涛 动 ;4M0O 为 大 西洋 
年 际 振荡 指数 


(a) SOI 和 SPEI 
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(c) SOI 和 PC2 


周期 /a 
Co 


1970 1980 1990 2000 2010 
年 份 


关系 数 均 通 过 显著 性 水 平 5% 的 检验 ,而 与 其 他 气 
候 因子 相关 系数 均 未 通过 检验 ,说 明 渭河 流域 干旱 
分 布 受 ENSO 影响 最 大 。 同 时 ,SPEI 与 SOI 的 相关 
性 最 好 ,PC1 次 之 , 且 随 着 模 态 的 增加 , PCs 与 SOI 
的 相关 系数 逐渐 减 小 ,相关 性 也 逐渐 降低 。 

进一步 ,分 别 进行 SOI 与 SPEI PCs 的 交叉 小 
波 变换 分 析 ,研究 时 间 序 列 彼此 之 间 的 时 频 变 换 关 
系 ,结果 如 图 7 所 示 。 图 中 被 粗 黑 线 包 围 区 域 表示 
通过 5% 显著 性 水 平 的 标准 红 噪 声 检验 ,箭头 方向 
反映 了 气候 因子 与 SPELT, PCs 的 相位 差 ,其 中 由 左 
向 右 的 箭头 表示 气候 因子 与 SPEL PCs 同 相位 ,由 
从 右 至 左 箭头 表示 反 相 位 ,而 垂直 向 上 或 向 下 的 箭 
头 分 别 表 示 气 候 因 子 超前 或 落后 SPEI PCs 时 间 序 
列 1/4 PAW, 

从 图 7a 和 7b 中 可 以 看 出 ,南方 涛 动 指数 SOT 
与 SPEI .PC1 与 交叉 小 波 变换 结果 相近 ,同样 可 以 
说 明 渭 河流 域 干旱 变化 主要 受 模 态 1 的 控制 影响 。 
SOI 对 PCL 的 影响 近年 来 主要 集中 在 9 ~ 13 a 的 区 
间 , 其 中 1995 一 2010 年 中 的 9 ~ 13 a 区 间 表 现 为 近 


(b) SOT 和 PCl 
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(d) SOT 和 PC3 
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图 7 SOI +3 SPEI PCs 的 交叉 小 波 变换 
Fig.7 Cross-wavelet transform between SOI and SPEI, PCs 


tó wy 


负 相 关 关 系 ,流域 的 干 湿 变 化 与 S01 指数 的 变化 相 
即 SOL 值 越 大 ,其 相对 应 的 SPEI 值 越 小 ,流域 遭 
全 流域 干旱 的 事件 越 多 (图 7b); SOL 对 PC2 的 
sd 4 ~7 a 的 区 间 , 其 中 1980—1995 
年 中 的 4~7 a 区 间 表 现 为 近 负 相关 关系 , 即 SOL 值 
越 大 ,流域 空间 结构 上 发 生 西 北 地 区 干旱 而 东南 地 
区 湿润 的 事件 将 会 越 多 (图 7c) ;S01 对 PCS 的 影响 
近年 来 主要 集中 在 8 ~ 14 a 的 区 间 , 其 中 1985 一 
2015 年 中 的 8 ~14 a 区 间 表 现 为 典型 的 负 相 关 , 即 
SOI 值 越 大 ,流域 空间 结构 上 发 生 东 部 地 区 干旱 而 
西部 地 区 湿润 的 事件 会 更 多 (图 7d) 。 


4 4 i 


基于 渭河 流域 气象 站 点 历史 实测 降水 量 和 气温 
数据 ,本 文 以 SPEI 作为 干旱 评估 指标 ,采用 M-K 趋 
势 检验 经验 正 交 函数 以 及 小 波 变换 等 方法 分 析 了 
渭河 流域 干旱 时 空 变换 特征 ,并 探讨 了 流域 干旱 与 
大 尺度 气候 因子 之 间 的 响应 关系 。 通 过 本 文 研究 得 
出 的 主要 结论 : 

(1) 在 过 去 1961 一 2017 年 共 57 a 间 ,渭河 流域 
整体 上 呈现 变 旱 的 趋势 。 尤 其 ,在 1990 年 之 后 流域 
干旱 程度 和 发 生 现象 明显 增加 , SPEL 平均 值 从 
1961 一 1990 年 的 0. 34 降低 到 1991 一 2017 年 的 
-0.37, 

(2) 通过 经 验 正 交 函数 分 析 表 明 ,渭河 流域 干 
湿 分 布 场 主要 有 3 种 典型 模 态 类 型 ,分 别 为 :全 局 分 
布 型 ,表现 为 渭河 流域 都 干旱 或 者 都 湿润 的 特征 ,对 


应 模 态 特征 方差 贡献 率 为 50. 43% ,是 渭河 流域 干 

温 分 布 的 最 主要 控制 模 态 ;西北 一 东南 反问 分 布 型 ， 
表现 为 流域 西北 地 区 干旱 东南 地 区 湿润 ,或 者 西北 
地 区 湿润 ,而 东南 地 区 干旱 ,该 模 态 对 应 特征 方差 贡 
献 率 为 17.35% ; 东 一 西 反 向 分 布 型 ,表现 为 东部 地 
区 湿润 西部 地 区 干旱 ,或 者 东部 地 区 干旱 西部 地 区 
湿润 , 模 态 对 应 特征 方差 贡献 率 为 5.06% 。 

(3) 在 6 种 常用 的 大 尺度 气候 因子 SOI, PDO, 

NP NAO AO 及 AMO 中 ,表征 ENSO 事件 的 SOI 对 
渭河 流域 干旱 时 空 分 布 影响 最 大 。 同 时 ,从 相位 上 
分 析 看 ,总 体 上 SOT 值 越 大 ,流域 发 生 干 旱 事件 相对 
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Abstract; The Weihe River Basin is located in the semiarid and semihumid transition region of northern China 
and experiences frequent drought events; therefore , study of droughts in the Weihe River Basin is of great impor- 
tance. Based on the precipitation and temperature data from 28 meteorological stations ,the spatial and temporal vari- 
ation characteristics of dryness were analyzed by using the Mann-Kendall ( M-K ) trend test, empirical orthogonal 
functions ,and wavelet transform analysis. At the same time, the relationship between drought characteristics and six 
large-scale climatic factors was studied, and the potential impact of the main factors on the spatial and temporal 
characteristics of drought was further discussed. According to the research, we obtained results showing that the 
Weihe River Basin exhibits a drought trend between 1961 and 2017. Particularly , after 1990 ,the drought degree and 
number of drought phenomena in the basin increased greatly and the average Standardized Precipitation Evapotrans- 
piration Index ( SPEI) decreased from 0.34 during 1961 - 1990 to -0.37 during 1991 - 2017. Empirical orthogo- 
nal function analysis shows that there are three typical modes in the wet and dry distribution field of the Weihe Riv- 
er Basin. The global type with corresponding modal eigenvalues of 50.43% is characterized by either drought or hu- 
midity and is the main control mode of dry distributions in the Weihe River Basin; northwest-southeast reverse dis- 
tribution type with corresponding modal eigenvalues of 17.35% shows that the northwest and southeast regions have 
different drought characteristics; east-west reverse distribution type with corresponding modal eigenvalues of 5.06% 
shows that the eastern and western regions have different drought characteristics. Finally , the research reveals that 
the Southern Oscillation Index (SOI) , which represents ENSO events , exhibits the strongest correlation with drought 
distribution in the Weihe River Basin and from the perspective of phase analysis ,abnormal SOI variations will gen- 
erally bring more drought events to the Weihe River Basin. 
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